Head-on/Near Head-on Collisions of Neutron Stars With a Realistic EOS 
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It has been conjectured that in head-on collisions of neutron stars (NSs), the merged object would 
not collapse promptly even if the total mass is higher than the maximum stable mass of a cold NS. 
In this paper, we show that the reverse is true: even if the total mass is less than the maximum 
stable mass, the merged object can collapse promptly. We demonstrate this for the case of NSs 
with a realistic equation of state (the Lattimer-Swesty EOS) in head-on and near head-on collisions. 
We propose a "Prompt Collapse Conjecture" for a generic NS EOS for head on and near head-on 
collisions. 

PACS numbers: 04.25.Dm,04.30.+x,97.60. Jd,97.60.Lf 



Introduction. Shapiro |l| has conjectured that two 
neutron stars (NSs) falling from infinity and collid- 
ing head-on would not, independent of their masses, 
promptly collapse to a black hole. The basic argument for 
the conjecture is that, until there is significant neutrino 
cooling, the thermal pressure generated by shock heating 
is always enough to support the merged object. Since 
neutrino cooling operates on time scales of seconds in- 
stead of the NS dynamical time scale of milliseconds, this 
could have significant implications on gravitational wave 
and neutrino emissions. The results obtained in 0] were 
based on a polytropic equation of state (EOS) P = Kp T 
(with K a function of the entropy), and it was suggested 
that the result was true for a general EOS. 

We have shown in Q that one implicit assumption 
used in the derivation in 0, namely that the collision 
can be approximated by a quasi-equilibrium process, is 
not valid. We carried out simulations of the head-on col- 
lision of neutron stars described by a polytropic EOS, as 
in the conjecture. For two 1.4 M Q NSs (with a poly- 
tropic index of T = 2, initial polytropic coefficient K of 
1.16 x 10 5 cm 5 /g s 2 as in a typical NS model), we showed 
that the merged object collapsed promptly. The shock 
front generated in the collision does not even have time to 
propagate to the outer part of the merged object before 
it is engulfed in an apparent horizon, let alone produce 
enough thermal pressure to support the merged object as 
envisioned in the quasi-equilibrium argument of [l|. In 
our study in Q{, the merged object had a mass well above 
the critical mass (the maximum mass that the EOS with 
the polytropic constants T and K given above can sup- 
port). We pointed out in 0,0 that the prompt collapse 
is due to the dynamical compression of the collision that 
is absent in the quasi-equilibrium analysis in Q, 0] . 

This brings up three further questions: 1. What if 
the mass of the merged object is less than the critical 



mass? Will the dynamical compression in the collision 
process be strong enough to initiate a collapse? 2. What 
if we use a realistic EOS instead of a polytropic EOS? 
3. What if we break the exact axisymmetry? In this 
paper we answer these three questions with one set of 
numerical simulations. The simulations are based on the 
GR-Astro code (formerly called GR-3D) constructed in 
the NASA Neutron Star Grand Challenge project [5j and 
the NSF Astrophysics Simulation Collaboratory (ASC) 
project 0. For the construction of the code and the 
classes of validation tests we have carried out for it, see 
The GR-Astro code solves the coupled set 
of the Einstein equations and the general relativistic hy- 
drodynamic equations with a realistic EOS. It will be 
released to the community through the ASC portal 
upon completion of the project. 

In this paper, we report on simulations of NSs con- 
structed with the Lattimer-Swestv|llj (LS) EOS, which 
has been used in various neutron star studies (all existing 
simulations based on the LS EOS that we are aware of 
arc based on Newtonian gravity and hence cannot answer 
questions of collapse). To the best of our knowledge, our 
simulations represent the first set of general relativistic 
3D simulations based on a realistic EOS (for polytropic 
EOS simulations, see 0). In this study, we use neu- 
tron stars of rest (baryonic) mass 1.6 Mq (corresponding 
to an ADM mass of 1.4 Mq in isolation), with a radius 
of 13.8 km (proper distance). The merged object has a 
rest mass of 3.2 Mq which is considerably lower than the 
critical mass of 3.67 Mq in the LS EOS with our choice 
of parameters. Nevertheless, we find that the merged 
object will promptly collapse to a black hole within a 
dynamical timescale. An apparent horizon is found en- 
gulfing the shock wave at 0.15 ms after the two stars have 
touched, with time measured at infinity, i.e. at the edge 
of the computational grid, which is 47 km away from the 
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collision center. The prompt collapse is verified for both 
the head-on collision case and an off-axis collision (the 
axisymmetry is broken by an impact parameter of 1/2 
stellar radius). This study demonstrates that dynami- 
cal effects are strong enough to cause a prompt collapse 
even when the total mass is below the single-star critical 
mass, and that this result is not a consequence of ex- 
act axisymmetry. We hence post the following "Prompt 
Collapse Conjecture" : For head-on and near head-on col- 
lisions of neutron stars described by a generic equation 
of state and infalling from rest at infinity, there exists 
a window in the rest mass of the merged object, below 
the critical single-star rest mass, where prompt collapse 
to a black hole can occur. The claim of the "near head- 
on" part means that the prompt collapse is stable with 
respect to small perturbations of the initial velocity. 

The Setup. In this paper we use our implementa- 
tion of the LS EOS as described in 0. We use the LS 
EOS in a tabular form with rest mass density and specific 
energy density as the two independent thermodynamic 
variables, which are evolved using the general relativis- 
tic hydrodynamic (GR-Hydro) equations (the lepton to 
baryon ratio is set to a constant (0.1) in all simulations 
in this paper). As in the Newtonian LS EOS simulations 
of [3 we set the initial specific energy density of 
the NSs at a relatively large value of 0.9 Mev. The ADM 
mass and rest mass as a function of central rest mass 
density for a single static NS is given in Fig. 1. We see 
that the critical rest mass of an LS EOS star is at 3.67 
Mq. In our simulation we use NSs with a rest mass of 
1.6 Mq (marked by a X in Fig. 1). The merged object is 
hence guaranteed to have a rest mass below the critical 
single-star mass. In Fig. 1 we have also plotted the cor- 
responding curves for the case of a polytropic EOS with 
r = 2 and K = 1.16 x 10 5 cm 5 /g s 2 for comparison. 

In the head-on collision case, we put the two stars at 
a proper distance of d = 42 km apart (about 3R separa- 
tion, R — radius of star) along the z-axis, and boost them 
toward one another at the speed (as measured at infinity) 
of y/GM/d (the Newtonian infall velocity). The metric 
and extrinsic curvature of the two boosted stars are ob- 
tained by (i) adding the off-diagonal components of the 
metric, (ii) adding the diagonal components of the met- 
ric and subtracting 1 , (iii) adding the components of the 
extrinsic curvature. The resulting matter distribution, 
momentum distributions, conformal part of the metric, 
and transverse traceless part of the extrinsic curvature 
are used as input to York's procedure [lfij ] for determin- 
ing the initial data. The initial data satisfies the com- 
plete set of Hamiltonian and momentum constraints to 
high accuracy (terms in the constraints cancel to 10~ 6 ), 
and represents two NSs in head-on collision. 

The initial data is then numerically evolved by solv- 
ing the coupled Einstein GR-Hydro evolution equations 
with numerical methods described in . The simulations 
reported here use the "1 + log" slicing ^t|- The simula- 
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FIG. 1: The ADM mass and rest mass vs. central rest 
mass density of a static NS with the LS EOS. The corre- 
sponding curves for a polytropic EOS with Y = 2, K = 
1.16 x 10 5 cm 5 /g s 2 are plotted for comparison. 



tions have been carried out with resolutions ranging from 
Ax = 0.74 km to 0.3 km (28 to 70 grid points across 
each NS) for convergence and accuracy analysis. We find 
that the constraint violations rise linearly throughout the 
evolution, and converge to zero with increased numerical 
resolution. The total rest mass of the system is conserved 
to better than 0.2% throughout the simulations. 

The Results. In Fig. 2a, we show the collapse of 
the lapse along the z axis from t = ms to t — 0.37 ms 
at intervals of 0.0926 ms. (With the reflection symme- 
try across the z = plane and the axisymmetry of the 
head-on collision, we only need to evolve the first oc- 
tant.) By the time t — 0.37 ms, the lapse has collapsed 
significantly. Fig. 2b shows the corresponding evolution 
of the zz component of the metric function. The "grid 
stretching" peak, characteristic of a black hole evolved in 
a singularity avoiding slicing, is apparent. 

Fig. 3 shows contour lines in the y — plane 
of the log of the gradient of the rest mass density 
log f ^/V*(/o)Vj(/ci)^ at time t — 0.37 ms. Sharp changes 
in rest mass density (where contour lines bunch up) in- 
dicate shocks. We see that the shock front is at 7 km 
in the x direction, and 10 km along the z direction (the 
collision axis) and has not yet reached the back end of 
star (at 14.5 km). At this point the shock front is still 
moving outward in coordinate space, although it is com- 
pletely engulfed by the apparent horizon, as seen in Fig. 4 
below. 

In Fig. 4, we show the intersection of the apparent hori- 
zon (AH) with the x — z plane. To confirm the location 
of the AH, convergence tests, both in terms of resolution 
and in terms of location of the computational bound- 



3 




0.4- 



0.2; 



— 1 = 


ms 


t = 


0.0926 ms 


.... t = 


0.185 ms 


— 1 = 


0.278 ms 


t = 


0.370 ms 




X 




1 5 



, km 



FIG. 3: Equally spaced contour lines of the log of the gra- 
dient of the rest mass density log ^ y^V^p)Vi(p)^ , showing 



the shock front at t = 0.37 ms. 
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FIG. 2: The lapse a and metric component g zz is displayed 
along the 2-axis at various times for the head-on collision of 2 
NSs using the LS EOS. This simulation used 163 3 grid points, 
with Ax = 0.3 km. 



ary, have been carried out. The solid, short dashed, long 
dashed and dotted lines correspond to the AH locations 
at resolutions of Ax = 0.74, 0.52, 0.37 and 0.30 km. The 
inset is a close-up view near the z-axis. We see that the 
location of the AH converges in a first order manner. We 
emphasize the importance of demonstrating the conver- 
gence of the final results in a numerical simulation. Such 
a validation of the final result is computationally expen- 
sive but nevertheless indispensable, as a short term con- 
vergence test of the numerical code does not guarantee 
the convergence of the final result at late time. In Fig. 5 
we show the hamiltonian constraint converging linearly 
with respect to resolution for the whole duration of the 
numerical evolution. 

Analysis of the Results. The above results demon- 
strate the phenomenon of a "sub-critical mass collapse" : 
a prompt collapse of the merged object, despite the fact 
that its total rest mass is below the critical rest mass of 
a single NS. 

One might question whether this result depends sen- 



FIG. 4: The position of the AH in the x — z plane at different 
resolutions, all at t = 0.37 ms. 



sitively on the choice of the initial velocity. As the star 
was only at about 3i? separation at the start of the sim- 
ulation the infall velocity determined by the Newtonian 
formula is not accurate. We have verified that increasing 
or decreasing the initial velocity by 10% does not change 
the qualitative feature of the collision, and in particular 
the prompt collapse result. Likewise, our initial configu- 
ration while satisfying the constraints may not have ex- 
actly the same gravitational wave content and distortion 
of the NSs as one that represents free fall from infinity. 
We do not expect that to have a significant influence on 
the result of prompt collapse (although it will affect the 
gravitational waveforms). To verify this point we have 
also carried out simulations using slightly distorted stars 
and slightly modified extrinsic curvature for the construc- 
tion of the initial data. We confirmed that the prompt 
collapse result is not affected. 

One might also wonder whether this prompt collapse 
result is a consequence of the exact symmetry used. To 
determine this point we carried out simulations with the 
two initial NS moving in the z direction but offset in 
the x direction by half a radius (an impact parameter 
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FIG. 5: The L2 norm of the hamiltonian constraint converg- 
ing linearly with respect to resolution for the whole duration 
of the numerical evolution. 
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FIG. 6: The position of the AH at t = 0.37 ms for a near 
head-on collision using a spatial resolution of Ax = 0.37 km. 



of 6.9 km). With this setup, the computation expense 
is significantly higher, as we can no longer evolve just 
an octant as in the axisymmetric case above. We found 
that the inclusion of a small impact parameter does not 
change the qualitative features of the collision, including 
the occurrence of prompt collapse. In Fig. 6, we show 
the apparent horizon found again at t = 0.37 ms. We see 
that the AH is tilted with respect to the z axis, but is 
very similar in shape and size to the head-on case. 

In all the above cases, we have prompt collapses. We 
have confirmed that at a low enough mass, the merged 
object will not collapse but will instead merge, bounce, 
oscillate and form a stable NS. However, the determina- 
tion of the dividing line between a final state black hole 
or a final state NS would require extremely high reso- 
lution whose computational expense is beyond what is 
available to our group. 

Finally, one may also ask whether the sub-critical mass 
collapse occurs in the head-on collision of NSs described 
by a polytropic EOS. Preliminary studies with our 3D 
code show that, in fact, the dividing line in the rest mass 



of the final merged object between prompt collapse and 
non-prompt collapse is quite close to the critical rest mass 
of a single NS for the T — 2 polytrope. We are currently 
working on an axisymmetric version of the code which 
we will use to not only explore this question, but to also 
explore the possible existence of type-1 critical phenom- 
ena at the interface dividing the prompt collapse and 
non-prompt collapse cases. 

Conclusions. In this paper we report on a sub- 
critical mass collapse phenomenon in the head-on/near 
head-on collisions of NSs with a realistic EOS. We pro- 
pose a "Prompt Collapse Conjecture" : For head-on and 
near-head-on collisions of neutron stars described by a 
generic equation of state and infalling from rest at in- 
finity, there exists a window in the rest mass of the 
merged object, below the critical single-star rest mass, 
where prompt collapse to a black hole can occur. This 
is the opposite of the Shapiro conjecture 0, which pre- 
dicted no prompt collapse for all NS masses, including 
those above the single-star critical mass. 

It has been argued in 01 that head-on/near head-on 
collisions of NSs could have a significant event rate, and 
could be a candidate for a sub-class of short gamma-ray 
bursts. The results of prompt collapse reported in this 
paper could have implications on the observation of such 
processes, with the prompt formation of the horizon cut- 
ting off the causal connection of the shock heated matter 
from outside observers. Note that we are not claiming a 
prompt collapse in inspiral coalescence of NSs. 
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